Quantum criticality in a generalized Dicke model 
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We employ a generalized Dicke model to study theoretically the quantum criticality of an ex- 
tended two-level atomic ensemble interacting with a single-mode quantized light field. Effective 
Hamiltonians are derived and digonalized to investigate numerically their eigenfrequencies for dif- 
ferent quantum phases in the system. Based on the analysis of the eigenfrequencies, an intriguing 
quantum phase transition from a normal phase to a super-radiant phase is revealed clearly, which 
is quite different from that observed with a standard Dicke model. 



INTRODUCTION 



Quantum phase transition (QPT) and quantum criti- 
cal phenomena, which are induced by the change of pa- 
rameters and are accompanied by a dramatic change of 
physical properties, occur at zero temperature in many- 
body quantum systems [jj]. Usually, a QPT may emerge 
in the parameter region where there is the energy level 
crossing or the symmetry-breaking. QPTs have been 
mainly studied in connection with correlated electron and 
spin systems in condensed matter physics [l|. Very re- 
cently, it has also been paid much attention in the light- 
atoms interacting systems |2|, y| , which enables us to un- 
derstand the transition from radiation to super-radiation 
from a different viewpoint. 

The systems of atomic ensembles interacting with opti- 
cal fields have been studied both experimentally and the- 
oretically, e.g., the electromagnetic induced transparency 
Q and the quantum storage of photon states [3, Q • The 
thermal phase transition phenomena [TJ have been stud- 
ied in the Dicke model 8] (that is, a two-level atomic en- 
semble coupling with optical field) or generalized Dicke 
models [9|. In particular, the QPT in a radiation-matter 
interacting system was recently explored based on the 
Dicke model, but merely with the single-mode Dicke 
model [2, [sj, |lO|, ll]|. When the coupling parameter A 
varies from that less than the critical value Ac to that 
larger than A^, the system goes from the normal phase to 
the super-radiant one in the presence of the symmetry- 
breaking. As the precursors of the QPT, the onset of 
chaos [2] and the entanglement properties [IQ] were stud- 
ied in detail. However, it is noticed that these stud- 
ies focused only on atomic ensembles with small dimen- 
sions compared with the optical wavelen gth , in which the 
dipole approximation can be used [2, y, llO| ■ In this spe- 
cial case, the light-atoms interaction is irrelevant to the 
spatial positions of atoms. But, generally speaking, a 
realistic atomic ensemble may extend in a large scale so 
that light-atoms interaction is spatially dependent [a.ll3|. 

In this paper, an exotic QPT phenomenon is investi- 



gated theoretically by developing the Dicke model for a 
more general case beyond the dipole approximation. We 
find that a kind of quantum critical phenomenon also 
occurs in this extended atomic ensemble, in which each 
atom interacts with a single mode quantized light field; 
but the quantum criticality is quite different from that 
deduced from the spatially independent Dicke model [2] ■ 
In the present study, a normal phase and four possible 
super-radiant phases are found, with only one of the four 
exhibiting the same critical point as that in the normal 
phase. Remarkably, it is shown that the ground-state 
energy in the above-mentioned superradiant phase, con- 
nects continuously to the normal phase one at the critical 
point, but its second drivative does not. 



A GENERALIZED DICKE MODEL FOR AN 
EXTENDED ATOMIC ENSEMBLE 

Let us consider an extended ensemble with N identical 
two-level atoms interacting with a single-mode quantized 
light field. Here, the spatial dimension of the atomic en- 
semble is much larger than the optical wavelength of the 
field. This radiation-matter system is usually described 
by a generalized Dicke model [8] with the Hamiltonian 
H = Hf) + Hi, 



N 



Ho = uja^a + ojo N^ cr. 



CC 7 



i=i 



N 



Hr 



Y^^e-^kr,^^^kr,^^in+^if))_ (1) 



i=i 



Here, h — 1, fc is the wave vector of the quantized light 

field and Vj is the position of jth atom; ace is the pop- 

(1) 
ulation operator of the j atom; (Teg is the flip operators 

between the excited state |e) and ground state \g) of the 

J atom with the same energy differences uq for all the 

atoms; a (a^) is the annihilation (creation) operator of 

the quantized light field, with A the relative coupling pa- 



rameter. For simplicity, the ensemble is assumed to be 
one-dimensional with its direction along the wave vector. 
Different from a standard Dicke model for small- 
dimension atomic ensembles y, the spatial-dependent 
factors exp(±ifcrj) are taken into account seriously 
though the momentum of the center of mass can be ne- 
glected. It is also remarked that the terms connected to 
the nonrotating-wave scenario are still kept in Hamilto- 
nian ([T]); in fact, if the rotating- wave approximation were 
used, the factors exp(±ifcrj) would be absorbed into (Teg 
and a 



M("), A("), B(") are defined as 
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NORMAL PHASE 

We can first introduce the following collective opera- 
tors %: 



ST 
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crli^e**'''^-. 
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It is obvious that in the limit of large A'^ with a small 
number of excitations (referred to as the normal phase) , 
namely, the excitation numbers in states |e) are much less 
than N ^ the above two operators approximately satisfy 
the independent bosonic commutation relations 



[i?,i?t]«[c,Ct]«l, [i?,Ct]«0 



(3) 



in the present extended ensemble, and can approximately 
be re-expressed as the two independent Bose operators 6^ 
and c^ 



[&,5t] = [c,ct]=l, [fe,ct]=0. 



(4) 



According to Ref. [l5|, we diagonalize the Hamiltonian 
dl]) in two steps: (i) find a unitary canonical transforma- 
tion T(") such that r(")^r(")tf^ ^ 1^ 7^(n)* ^ 7T(")7, 
where 
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^=1 -1 ' ^- VI 
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and (ii) introduce the quasiparticle operators 

V(») = r(")u(") - (/ii, h^, ha; hi hi hlf. 
Then the Hamiltonian ^ is cast into a diagonalized form 

H(-^=J24'\hlh. + \)-\trA(-\ (9) 

and describes the quasiparticle excitations with frequen- 
cies 1^12 3' which are obtained by diagonalizing fjM''^^^ 
with r(") into 



Tf]MT-^ ^ n = 



-Lu r 




(10) 



In the present normal phase case, the original radiation- 
matter system described by Eq. ^ is approximated as 
a coupling three-mode bosonic system with the "low en- 
ergy" effective Hamiltonian 



+Aa^(ct +6)+H.c.. 



(5) 



We now apply a Bogoliubov transformation to diago- 
nalize the above quadratic Hamiltonian ^ [l5|. First, 
we rewrite it as 

^(„) ^ iu(«)tM(")u(") - itr A("), (6) 

2 2 ^ ^ 

where the operator-valued vectors U*-"^ and the matrices 
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FIG. 1: (Color online) The real part (thin solid lines) and 
imaginary part (thick dashed lines) of the eigenfrequencies vs 
the coupling parameter A in units of ujo in the normal phase 
at the resonance case a; = cjq. 

Below, we focus only on the properties of eigen- 
frequencies in order to explore the existence of quantum 



criticality. Since the matrix M^"^ is of 6 x 6, general 
analytic results for the diagonalization of fjM^^' are dif- 
ficult to obtain. Nevertheless, we can diagonalize it nu- 
merically to obtain the eigen-frequencies ^-^ 2 3 UM- The 
related canonical transformation matrix T*-"-* can also be 
obtained numerically. Figure [1] shows the numerical re- 
suits for the real and imaginary parts of uji 2 3- For sim- 
plicity, we illustrate the resonant case to = luq — I here. 
Certainly, the non-resonance cases can also be studied 
numerically, with similar features being revealed. As seen 
from Fig. [2 when A > 0.7698, the imaginary parts of two 
eigen-frequencies are non-zero. This means that the cor- 
responding eigen-frequencies are complex and thus the 
eigen state is unstable and physically impossible. But it 
is inappropriate to consider naively that A — 0.7698 as 
the critical point. Since the eigenvalue Wg is negative 
in the range of A € (0.7071,0.7698), a negative eigen- 
frequency of the boson-mode is not allowed physically 
either. Therefore, in the resonance case uj — luq {— 1), a 
real critical point is located at A^"^ = 0.7071 = \/2/2. In 
addition, for a general case, the critical point is found to 
be A(") = ^/ZJS^. 



SUPER-RADIANT PHASE 

In order to describe excitations in the parameter re- 
gion above the critical point, we now incorporate the fact 
that both the field and the atomic collective excitations 
acquire macroscopic occupations, namely, the above ap- 
proximation to neglect the number of excitations over 
iV is no longer valid [2|- To this aspect, the introduced 
collective operators H^ and C"^ in Eq. ^ should be ex- 
pressed approximately as [l7| 



5^-^Vi-^, 



(11) 



^^-Vi-f' 



in terms of the Bose operators &^, 6, c^ and c. This trans- 
formation maps the original light-atoms system to a cou- 
pling three-mode bosonic system with the Hamiltonian 



H 



{^) - 



= uJo{wb + c'c) + uja'a + Xa* 



(12) 




H.c. 



where a, (3 and 7 are generally complex parameters in the 
order of 0{\N) [4I to be determined later. This is equiv- 
alent to assume that all modes behave as the nonzero, 
macroscopic mean fields above AJ;"^ . 

Keeping the terms up to the order of 0{N^), the 
Hamiltonian p^ becomes 



hm = 




■i'/'+ '°'^'„:^'-"' h/'+7-/) 



[ae^ 



2fcf 

-p*d) 



Zfee 



(/3et + [3*e) 



h.c. 



+Lud^d + LOec'^e + LOff'^f + Cq 



(s) 



(14) 



where the constant term 



-,('*) 



uj\a\ -fcJo( 
_2Aa 



f + l7l') 



fco + 7V ^f ) 



(15) 



will substantially contribute to the ground state energy 
at critical point; the renormalized frequencies 



uje = luq + Xa* j3/\J Nkc, 



UJO + Xaj/y/Nke 



with fco = N— |/3| and k[ — N— {-fl" . In the derivation of 
Eq. p4p , the terms being linear in the bosonic operators 
are eliminated by choosing the appropriate displacements 
a, (3 and 7 in the following four cases: 



v(1.2,3,4) _ 



pi3) ^ e^0 



M) 



or 




(16) 



X^, (j = l,3), 
/3(^' = e'^yxl, (j=2,4). 
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(i) 



2,4), 



where X± 



N 



2 (1 ± t^), and 



is an arbitrary real 
number relating to the phases of displacements. In fact, 
we see from the form of a^^' in Eq. (J16p that only when 



4A^ 



lu'^ujq > 0, 



(17) 



a^^' can be physically meaningful. Thus in the following 
discussions, it is required that 



A^ 



^ (. aH). 



For the present super-radiant phase, the bosonic modes 
may be displaced in the following way: 



(13) 



at - 


^ d^+a*; 


6t - 


-^ et-/3*; 
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- P~l*, 



It is interesting to note that this threshed is just the 
critical point determined in the normal phase case. 

Since ij'''^((/)) in Eq. p4| can be transferred to a (j)- 
independent Hamiltonian iJ*-*' (0 = 0) through a unitary 
transformation 



UW 



J4,[dU+e-^e-f'<f) 



we need only to look into the spectra of i?(*)(0) in the 
four cases specified by Eq. p6|) . respectively. Because 
H'-'^^O) is quadratic in each case, which is diagonalized 
by using the same method presented above as 



3 



ijO) = 2^uj 



U) 



jmj^) 



itr A(^) + c(^\ 



(18) 



for the four cases j = 1,...,4. Here, the quasi-particle 
excitation is described by the boson vector operators 

e — (til ,^2,^3,^]^ , ^2 , ^3 ) 

in the j-th phase, where 



S^^ =a- 



vO) 



M) =b + 



U) 



, /(^") = c 



v(j) 



according to Eq. (fT3| . T^^' is still the introduced unitary 
transformation to diagonalize 



where 



^(1.3) ^ 
5(1.3) 






W A^ C^: 
















B^ 




with 



= 2X^X±/Nu;, 

= X{X^ - X±/2)/,/NX^, 

= -XX±I2^NX^. 



UJ± 

B± 

A± 
C± 

Clearly, ujI''' {i — 1,2,3) is the i-th eigenfrequency for 
the Hamiltonian i?^-'^ . Note that the canonical transfor- 
mation matrix T'^^ can be obtained numerically in the 
numerical diagonalization of fjM^^' . 

The numerical results of eigen-frequencies uji2 3 '^''• 
the coupling parameter A are plotted in Fig. [21 The 



curves for both the real and imaginary parts of eigen- 
frequencies of H^^' in the resonant case are shown Fig. 
[^a). It is found that the eigenfrequencies are physically 
reasonable when A > 0.7071 since the imaginary parts 
of all the eigenfrequencies are zero. This means a novel 
"quantum phase" emerges above the critical point 



Al^) 



UJUJQ 



0.7071. 



It is seen from Fig. [SJa) that the eigenfrequency cjg 
is always zero above A^"'^ , which implies that H^^^^ is re- 
duced to have two independent boson modes. It is re- 
markable that the critical point is just the same one as 
that determined in the normal phase AJ; ' — X^J^' — Ac, 
demonstrating the consistency of our analysis. In Fig. 




FIG. 2: (Color online) The real part (solid lines) and imag- 
inary part (dashed lines) of the eigenfrequencies vs the cou- 
pling parameter A in units of ujq in the super-radiant phase 
at the resonance case u = ujo, where (a,b,c,d) correspond 
respectively to the cases (1,2,3,4) specified in the text. 

[2jb) [Fig. El^d)], the numerical results for the real and 
imaginary parts of three eigenfrequencies of H^^^^ [H'-^^] 
in the resonant case. As seen from Fig. ^h) [Fig. [2Id)], 
only when 

A > A^^^ = 0.8112 (A^-*) = 0.8112), 

another possible "quantum phase" may appear as the 
imaginary parts of all the eigenfrequencies are zero. 
While for H^^\ as seen from Fig. Efc), only when 
A > 0.8457 = Aj, , the imaginary part of the all eigenfre- 
quencies are zero, indicating a possible "quantum phase." 

In Fig. [3ja), we plot together the eigenfrequencies vs A 
for the normal phase and first super-radiant phase. The 
eigenfrequencies oj" and uj\ (for i = 1, 2, 3) are contin- 
uous at the critical point, respectively. Comparing with 
the results in the spatially independent Dicke model [4|, 
our numerical studies show clearly that the excitation 
energy Wg in the normal phase vanishes as [A — Ad^" 

and the characteristic length scale I3 = l/yoj^ diverges 
as [A — Ac I at the quantum transition point Ac, with 
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FIG. 3: (Color online) (a) The eigenfrequencies for normal 
phase and the first super-radiant phase, (b) The ground state 
energy densities (A^ — fO^) for the normal phase and the 
super- radiant phases at the resonance case uj — luq — 1. 



the exponents given by w = 1/2, z — 2 on resonance; 
however, it is interesting to note that no critical expo- 



nents for u 
ified since 



,(1) 



,(1) - 



in the super-radiant phase can be spec- 
0. Meanwhile, for the ground state, 



(fltfl) /N — below Ac, while 



N 



v(l) 



/iV ex (A - Ac) 



above Ac, i.e., the field is macroscopically occupied. So 
a^^' may be understood as a kind of order parameter of 
the super-radiant phase, whose critical exponent is 1/2 
above Ac. In addition. Fig. 3(b) presents the ground- 
state energy densities as a function of coupling for all the 
possible phases. Clearly, the ground state energy densi- 
ties of the first super-radiant phase is always the lowest 
one above Ac, while the other three approach to it in the 
large A limit; moreover, it connects continuously with 
that of the normal phase but possesses a discontinuity in 
its second derivative at Ac through a detailed numerical 
analysis. From this viewpoint, together with the fact that 
the same critical point is determined from both sides of 
the normal phase and the first one, it is most likely that 
only the first super-radiant phase is a real physical one. 



Although the effect of non-RWA terms may normally 
be negligibly small, the present work (also see Ref. [2]) il- 
lustrates that it plays a meaningful role when the atomic 
number N is large, e.g., the criticality differs from that 
with the RWA. On the other hand, for actual atoms that 
may not be pure two-level ones, other atomic transitions 
may occur and spoil the present model before the non- 
RWA terms become important. Nevertheless, the present 
study is still theoretically interesting and valuable, par- 
ticularly relevant to some atomic systems (or artificial 
and atomic-like ones) wherein the energy spacing of any 
other transitions is much larger than that of the consid- 
ered two levels (or other transitions do not exist). For 
example, for a Dicke-like model consisting of many 1/2 
spins coupled to single mode bosonic field (by electrical 
dipole coupling- like type), other transitions do not exist 
in the spin systems. Then the counter rotating terms 
play an important role the when the coupling parameter 
is close to the critical value. 

In conclusion, based on a generalized Dicke model, 
we have investigated theoretically the quantum critical- 
ity of an extended atomic ensemble with a larger spa- 
cial dimension comparable to the optical wavelength of 
a quantized light field. A useful formalism is developed 
to study numerically eigenfrequencies of the system in 
different quantum phases. Comparing with the critical 
phenomenon around the critical point Ac — y/ujLOo/2 for 
atomic ensemble of small dimension 2], a rather differ- 
ent quantum criticality is revealed around the transition 
point (Ac — y/ujLUo/2 = v^Ac) from the normal phase to 
the super-radiant phase. 
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REMARKS AND CONCLUSIONS 

Before concluding this paper, we wish to remark briefly 
on the origin of the occurred QPT in the present work. 
From Fig. 3(a), it is clear that the energy level of the 



first excited state of the system (wg + Eq ) touches 



(«) 



?M^ 



(n) 



the ground state energy level Eq ' (or £'q ) at the crit- 
ical point. Obviously, it is this level touching that ac- 
counts for the emergence of the QPT and the correspond- 
ing quantum criticality in the present generalized Dicke 
model. It is also remarked that the A^ terms (where A 
is the vector potential) has been neglected here, as done 
in several pre vious works [3, y, |ll| , while the absence of 



A^ terms [18|, |19| seems to be crucial for the observed 
quantum phase transition in the present model, namely, 
the presence of A^ terms in the model Hamiltonian leads 
to vanishing of the criticality. 
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